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ABSTRACT 

We present a new formation scenario for NGC 2419, which is one of the most lumi- 
nous, one of the most distant, and as weh one of the most extended globular clusters 
of the Milky Way. We propose that NGC 2419 is the remnant of a merged star cluster 
complex, which was possibly formed during an interaction between a gas-rich galaxy 
and the Milky Way. To test this hypothesis, we performed numerical simulations of 27 
different models of star cluster complexes (CCs) moving on a highly eccentric orbit in 
the Galactic halo. We vary the CC mass, the CC size, and the initial distribution of 
star clusters in the CC to analyze the influence of these parameters on the resulting 
objects. 

In all cases, the vast majority of star clusters merged into a stable object. The 
derived parameters mass, absolute V-band magnitude, effective radius, velocity disper- 
sion and the surface brightness profile are, for a number of models, in good agreement 
with those observed for NGC 2419. Despite the large range of CC sizes, the effective 
radii of the merger objects are constrained to a relatively small interval. A turnover 
in the rgfi vs. Mend space leads to degenerate states, i.e. relatively compact CCs can 
produce an object with the same structural parameters as a more massive and larger 
CC. In consequence, a range of initial conditions can form a merger object comparable 
to NGC 2419 preventing us to pinpoint the exact parameters of the original CC, which 
formed NGC 2419. 

We conclude that NGC 2419 can be well explained by the merged cluster complex 
scenario. Some of the recently discovered stellar streams in the Galactic halo might 
be related to the parent galaxy, which produced the cluster complex in our scenario. 
Measurements of the proper motion of NGC 2419 are necessary to prove an association 
with one of the stellar streams. 

Subject headings: Galaxy: evolution — globular clusters: individual (NGC2419) — 
methods: numerical 
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INTRODUCTION 



Globular clusters (GCs) are among the oldest radiant objects in the universe. They have 
typically compact sizes with half-light radii of a few pc and absolute V-band magnitudes between 
My = —5 to My = — 10 corresponding roughly to masses between Mq and 10^ (see 



Brodie &: Stradeiil2006l . and references therein). 



The Milky Way has a rich GC system containing 150 GCs (|Harrislll996l ). Only 13 GCs (or 9%) 
have an effective radius larger than 10 pc (see Figure [1]). Most of these extended clusters (ECs) 
are fainter than about My = —7, only NGC 2419 has a high luminosity of about My = —9.4 mag. 
While all GCs with a comparable luminosity are quite compact, NGC 2419 has a projected half- 
light or effective radius of about 20 pc. The structural properties of NGC 2419 derived by various 
authors are displayed in Tabled) NGC 2419, locate d at (/,6)=(1 80?4, 25?2), is one of the most 
metal-poor GCs with a metallicity of [Fe/H] = -2.12 (lHarriElll996l) at a Galactocentric distance of 



about 91.5 kpc and a heliocentric distance of 82.4 kpc (IHarris Il996l ). NGC 2419 shows a single 



stellar population with an age of 12 to 13 Gyr (|Salaris fc WeissI 



200 



FigureHalso show s the GCs from the LM C dMackev fc Gilmorel[20o3 : Ivan den Bergh &: MackevI 



20041'). M51, and M81 (jChandar et al.ll2004l ). from 68 dwarf galaxies of nearby loose ga laxy groups 
(jGeorgiev et al.ll2009l). and the nearby elliptical galaxy NG C 5128 (IHarris et al.ll2002l). as well as 



the ECs from M31 (jHuxor et al.l l2004l : iMackev et al.l \2Q0(i). M33 (IStonkute 



et al 



20081), Scl-dEl 



(|Da Costa et al.l l2009h . and EC 90:12 from NGC 1399 (|Richtler et al.l bood ). It is evident that 
NGC 2419 has a rather isolated location in the My vs. reS space. 

The unusual size and mass of NGC 2419 led to the speculatio n that NGC 2419 is not a GC, bu t 
the stripped core of a dSph galaxy captured by the Milky Way (e.g. Ivan den Bergh &: Mackevll2004l ) . 



A stripped core should, howeve r, show multiple stellar populations from the progenitor galaxy. 
In contrast, iRipepi et al.l (120071 ) demonstrated that NGC 2419 shows no sign of m ultiple stellar 



copulations and concluded that NGC 2419 is not a stripped core of a dwarf galaxy. iRipepi et al 



(|2007j ) also found that NGC 2419 is an Oosterhoff II cluster. They concluded that an extra-galactic 
origin of NGC 2419 is unlikely as the GCs of the Galactic dwarf galaxies usually show parameters 
between Oosterhoff types I and II. On the other hand, differences between NGC 2419 and typical 
values of GCs associated with Galactic dwarf spheroidal galaxies do not rule out an extra-galactic 
origin of NGC 2419 per se. For instance, five out of twelve GCs from the LMC are Oosterhoff II 
clusters (jCatelan Ibooa ). 



High-resolution imaging of gas-rich galaxies experiencing major interactions using the Hubble- 
Space Telescope has resolved regions with very intense star formation bursts. The bursting regions 
are typically located within the severely perturbed disks or tidal tails and are constrained to small 
complexes that contain a few to hundred s of young massive star clusters. Examples of such sy stems 
are the knots in the Antennae galaxies (IWhitmore &: Schweizerl Il995l: IWhitmore et al.lll999l). th e 
comple xes in the NGC 7 673 star-burst (jHomeier et ahlbood ). M8 2 dKonstantopoulos et ahl bood). 
Arp24 (|Cao fc Wul[2007l ). the "bird's head galaxy" NGC 6745 dde Griis et al.ll2003l ). NGC 6946 
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(jLarsen et al.lbooa '). Stephan's Quintet (jGallagher et alJboOlh . and NGC922 (jPellerin et al.lboid ') 

The masses of such complexes vary from 10^ Mq up to a few 10® Mq 



Bastian et al.1 (|2006l ) 



observed star cluster comple xes (CCs) in the Ante nnae with masses of the order ~ 10^ Mq and sizes 
of the order 100 to 200 pc. iPellerin et al.1 (120101 ) found young massive CCs with masses between 
10*^ Mq and lO'^'^ M© and diameters between 600 pc and 1200 pc in the collisional ring galaxy 
NGC922. 



Mengel et al.1 (120081 ) observed individual young (~10 Myr) star clusters associated with cluster 
complexes in the Antennae and NGC 1487. They compared dynamical mass estimates with derived 
photometric masses and found them in excellent agreement, implying that most of them survived 
the gas removal phase and are bound stellar objects. These young clusters are sufficiently stable to 
be used as building blocks for numerical simulations. iBastian et al.l (120091 ) found three 200 to 500 
Myr old, apparently stable clusters in the Antennae with very high radial velocities relative to the 
galactic disk, indicating that these star clusters will most likely become future halo objects. One 
cluster is surrounded by so far unmerged stellar features in its vicinity. 

It has already been shown that CCs can merge to form a variety of spheroidal stellar-dynamical 
objects, such as ultra-compact dwarf galaxies (UC Ds ), faint fuzzies and possibly dwarf spheroidal 
galaxies dKroupal Il998l : IPellhauer fc Kroupal l2002al lbl: ] Briins et al.l l2009l). In particular, the yo ung 
UCD W3 is most naturally understood to be a merged massive CC (IFellhauer &: Kroupal l2005l ). 



In the present paper, we apply these ideas to NGC 2419 by analyzing how the structural 
parameters of the final merger objects in the outer Galactic halo correlate with the underlying CC 
parameters and compare them with the observed parameters of NGC 2419. 

In Section [21 we describe the method and the parameters used for the calculations. Section [3] 
presents the results of the simulations, which will be discussed in Section HI In Section [5] we provide 
a summary. 



2. NUMERICAL METHOD AND SET-UP 

The formation scenario described in this paper starts with newly born complexes of star clusters 
in the Galactic halo. We model the dynamical evolution of various CCs leading to merger objects. 
We do not, however, consider the galaxy-galaxy interaction, which formed the CCs in the first 
place. 



2.1. The Code Superbox-H- 

SuPER BOX++ developed by Metz ( 20081) is a new implementation of the particle-mesh code 
SuPERBOX (|Bien et al.lll99ll : IFellhauer et al.ll2000l ). Differences in the implementation of the algo- 
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rithm make the new code Superbox++ much more efficient than the former. SuPERBOX had been 
implemented in the FORTRAN language with a particular focus on the minimization of usage of 
the random access memory (RAM), which is no longer a big issue in current days. Superbox-H- 
is now being implemented in the modern C-H- programming language using object oriented pro- 
gramming techniques. The algorithm has been implemented with a focus on the performance of 
the code, but at the same time keeping memory consumption at a low level. SuPERBOX+f makes 
particular optimal use of modern multi-core processor technologies. 

The code solves the Poisson equation on a system of Cartesian grids. The local universe is 
covered by a fixed coarse grid which contains the orbit of the CC around the center of the Milky 
Way. In order to get good resolution of the star clusters two grids with high and medium resolution 
are focused on each star cluster following their trajectories. The individual high resolution grids 
cover an entire star cluster, whereas the medium resolution grids of every star cluster embed the 
whole initial CC. All grids contain 128^ grid cells. The CC orbits in an analytical Galactic potential 
(see Sect. 12. 2. ID . For each particle in the CC the acceleration from the galactic potential is added 
as an analytical formula to the grid-based acceleration computed by solving the Poisson equation. 
The coordinate system is chosen such that the disk of the Milky Way lies in the x-y-plane with 
origin at the Galactic center. 



2.2. Orbit of NGC 2419 



2.2.1. Gravitational Potential of the Milky Way 

In our computations the Milky Way is represented by an analytical potential, which consists 
of a disk-, a bulge-, and a halo component. 



The disk is modeled by a Miyamoto-Nagai potential (jMivamoto &: Nagailll975l ) 

^>disk(^,^;) = - 



(1) 



^2^2 



with Md = 1.0 X 10^^ M(7), Od = 6.5 kpc, and 6d = 0.26 kpc. The bulge is represented by a Hernquist 
potential (jHernquisti 1 19901 ) . 



bulge 



r + ab' 



with Mb = 3.4 X lO^'^ Mq, and Ob = 0.7 kpc. The halo is a logarithmic potential, 

\ 

with vq = 186.0 km s^^, and rhaio = 12.0 kpc. 

This set of parameters gives a realistic rotation curve for the Milky Way. 



(2) 



(3) 
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2.2.2. Orbital Parameters 



The calculation of an orbit for NGC 2419 requires, next to the spatial coordinates, a good 
knowledge of the actual velocity vector of the cluster. As no proper motion measurements for 
NGC 2419 are available, the orbit cannot be properly fixed without major assumptions. 



Kind (|l962l ) suggested that the tidal radius of a GC is determined by t he or bital position with 
the highest gravitational force, which is the perigalacticon. Ivan den Berghl (|l995l ) used this method 
to estimate a perigalactic distance of 20 kpc for NGC 2419. While this method provides only a very 
rough estimate of the perigalactic distance, we use it to narrow down the range of possible orbits. 



NGC 2419 is one of the outermost Galactic GCs at a distance of i?gai = 91.5 kpc (|Harrislll996l ). 
The ver y large distance and the observed low (heliocentric) negative radial velocity of Whci = ^20. 3 
km s~^ (jBaumgardt et al.ll2009l ). which corresponds to a Galactic-standard-of-rest velocity of fcsR 
= -26.8 km s~^, suggests that NGC 2419 is close to its apogalacticon. 



We choose a velocity vector for NGC 2419 in a way that it 



• is consistent with the observed radial velocity, 

• produces an orbit between Galactic radii 20 and about 92 kpc, 

• is neither a polar orbit nor an orbit within the Galactic plane. 

The velocity vector (ux, fy, v^) = (45 km s^"*^, -56 km s~^, 47 km s^"*^) will be used as the current 
velocity vector of NGC 2419 for calculating the orbit. 

Figure [2] shows the orbit of NGC 2419 calculated back in time using the Milky Way potential 
as given in Sect. 12.2.11 The orbital period is about 1.3 Gyr. The CC will most likely be formed 
during the first perigalactic passage of the parent galaxy. However, as the calculated objects will 
not change significantly after about 7 Gyr of evolution (see Sect. 13. ip . we do not start at t = -12 
Gyr (the age of NGC 2419) in order to cut down the computing time. For all numerical simulations 
we choose the perigalactic passage at t = -9.568 Gyr as a starting point. 



2.3. Initial Configuration of the Cluster Complex and Model Parameters 



The CC models consist o f N^ ^ 



distribution flPlummer 



1911 



20 star clusters and are diced according to a Plummer 



Kroupall2008l ). The cutoff radius, R^^, of the CC is four times the 



Plummer radius, -Rpj*^. The initial velocity distribution of the CC models is chosen such that the 
CC is in virial equilibrium. A detailed description of the ge neration of init i al coo rdinates (space 
and velocity) for Plummer models is given in the appendix of lAarseth et al.l (jl974l ). 



The individual star clusters building up the CCs in our simulations are Plummer spheres with 
a Plummer radius of = 4 pc and a cutoff radius of ii^iS = 20 pc. Each star cluster has a mass 
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of M^^ = 0.05 X M^^ and consists of TVq = 100 000 particles. The velocity distribution of the 
individual star clusters is chosen to be initially in virial equilibrium. 

In total, we considered 27 different models (see Tables [2] and [3]) , which are denoted M_x_y_z, 
where x is the number of the initial configuration, i.e. the detailed distribution of the individual 
star clusters in the CC, y is the CC mass, M'^'^, in units of 10^ M©, and z is the CC Plummer 
radius, R^i", in pc. Figure [3] visualizes the CC parameter range covered in the M*^^ vs. R^^" 
space. Figure U] illustrates the different initial distributions. Figure 0^ and b are the same initial 
distribution of star clusters that were scaled according to their R^i", while Figure shows a less 
concentrated distribution of star clusters. 



3. RESULTS 

We carried out 27 different numerical simulations to get an estimate of the influence of varying 
initial CC conditions. All calculations start at the perigalactic passage at io = -9.568 Gyr and are 
calculated up to the current position of NGC 2419. 



3.1. Time Evolution of the Merging Process 

The merging process of model M_l_1.5_100 is shown in Figure [5] as contour plots on the xy- 
plane to illustrate the detailed evolution of the merging process. The snapshots were taken at f 
= t - to = 0, 50, 100, 300, 760 and 1500 Myr. At t' = 50 Myr the merger object is already in 
the process of forming, but the majority of star clusters are still individual objects. In the course 
of time more and more star clusters are captured by the merger object. Thus the merger object 
becomes more extended. After 10 crossing times (f = 760 Myr) there are still 2 unmerged star 
clusters in the vicinity of the merger object. In the last snapshot at t' = 1500 Myr the merging 
process is completed and 19 out of 20 star cluster have merged forming a smooth extended object. 
One star cluster escaped the merging process. It follows the merger object on its orbit around the 
Milky Way at a distance of about 14 kpc (at f = 9.568 Gyr). 

The timescale of the merging process depends on the initial CC mass, the CC size and the 
distribution of star clusters within the CC. For model M_l_1.5_100, 50 percent of the clusters have 
merged after approximately 100 Myr. The time within which half of the clusters merge increases 
with the CC size (from 15 Myr for model M_l_1.5_25 to 200 Myr for model M_l_1.5_150) and 
decreases with CC mass (from 150 Myr for model M_l_1.0_100 to 65 Myr for model M_l_3.0_100). 

The further time evolution of the merger object of model M_l_1.5_100 is plotted in Figure [6j 
The enclosed mass of the merger objects is defined as the mass within 800 pc. The half-mass radius 
is the radius of the sphere within which half of the mass is enclosed. However, as observers can 
only derive a projected half-mass radius, we calculate also a value projected on the sky defined 
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as the projected radius within which half of the mass is included. The projected half-mass radius 
is slightly smaller than the three-dimensional half-mass radius (Table [3|) and corresponds to the 
observed effective radius, r^s, if mass follows light. The effective radius and mass of the merger 
object decrease and become fairly constant after about 7 Gyr. The structural parameters change 
only very slightly in the next few Gyr. 

3.2. Variation of initial CC mass and size 

We consider models with three CC masses, M^"^ = 1.0, 1.5, and 2.0 xlO^ Mq, and six CC 
Plummer radii between R^y" = 25 and 150 pc in steps of 25 pc to analyze the dependence of the 
structural parameters of the merger objects on the initial CC mass and size. In addition, two 
models with i?^c ^ 2OO and 300 pc are calculated for M^*-^ = 2.0 x 10^. All models have the same 
relative initial distribution of star clusters, which is scaled according to the respective Plummer 
radii (see Figure 3^ and b). The velocities of the individual star clusters are also scaled with CC 
mass and size to keep the CCs initially in virial equilibrium. 

For all models the merging process leads to a stable object. The number of merged star clusters 
is between 10 and 20. The properties of the merger objects at the current position of NGC 2419 
are displayed in Tabled 

Figure [Tti shows the enclosed mass. Mend, of the merger objects as a function of the initial 
CC Plummer radius, R^i^- The fraction of the initial CC mass, which is bound to the merger 
object, decreases almost linear with increasing CC size from about 92% at R^f^ = 25 pc for all 
three masses to values of 46, 53, and 59% at R^^ = 150 pc for CC masses of M^*-^ = 1.0, 1.5, and 
2.0 X 10^, respectively. Mass-loss occurs either by the escape of individual stars from the diffuse 
stellar component, which builds up during the merging process (see Sect. 13. ip . or by entire star 
clusters escaping the merging process. For all three masses, all 20 star clusters merge for compact 
models up to R^^" = 75 pc, i.e. mass-loss is only from the diffuse component for these models. For 
the more extended models with R^p > 100 pc star clusters escape from the merging process. For 
extended models with R^^" = 150 pc six, five, and two star clusters escape for CC masses of 
= 1.0, 1.5, and 2.0 x 10^, respectively. 

Figure [7)3 shows the effective radii, r^g, of the merger objects as a function of the initial CC 
Plummer radius, R^i"- The effective radii of the merger objects increase with increasing CC size 
up to R^Y^ = 75 pc for M'-^*-' = 1.0 x 10^ and up to R^^^ = 100 pc for the more massive models. 
For larger values of R^^" the effective radii decrease rapidly. Despite the large range of CC sizes 
(i?^i^ = 25 to 300 pc), the effective radii of the merger objects are constrained to a relatively small 
interval, e.g. models with M'~'*-' = 2.0 x 10^ are constrained to values between 13.2 and 25.6 pc. 

Figure [8] combines Figures [7^ and b showing the effective radius of the merger objects as a 
function of their enclosed mass. Only the six radii R'^^^ = 25, 50, 75, 100, 125, and 150 pc, where 
data for all three masses are available, were used to allow for an overview on the trends (models 
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with the same initial CC Plummer radius are connected by dotted hnes, models with the same 
initial mass with solid lines). For all CC Plummer radii the effective radii increase with the initial 
CC mass. While the effective radii vary only slightly with mass for models with = 25, 50 and 
75 pc, the slope is considerably steepening for larger CC sizes. 

The formation process of the merger objects depends on the compactness of the initial CC. A 
measure of how densely a C C is filled with star cl usters for an equal number Nq''^ of star clusters 



is given by the parameter a (jFellhauer et al.ll2002l ). 



« = ]|c' (4) 

where i?pp and R^f" are the Plummer radius of a single star cluster and the Plummer radius of the 
CC, respectively. In our models, the individual star clusters are not scaled with R^^- This leads 
to decreasing values of a from 0.16 for R^f^ = 25 pc to a = 0.013 for R^^^ = 300 pc. In general, 
large values of a accelerate the merging process because the star clusters already partly overlap in 
the center of the CC, whereas small values hamper the merging process. 

Also the tidal field counteracts the merging process. An estimate of the influence of the tidal 
field on the CC is given by the parameter 

r>CC 

/3 = -# (5) 



( Fellhauer et al.ll2002l ). which is the ratio of the cutoff radius R^^ of the CC and its tidal radius 
rp*-'. The tidal radius, r^'^, is defined as the radius where the attractive force of the CC on an 
object equals the pulling force of the Milky Way. Using the Milky Way potential as given in Section 
12.2.11 and the CC potential, which is given by the Plummer potential 

the tidal radius at the peri-galactic distance of 20 kpc is r^^^ = 296, 291, 283, 274, and 262 pc for 
M^*-^ = 1.5 X 10^ and R^^" = 50, 75, 100, 125 and 150 pc, respectively. The corresponding values 
of /3 are 0.68, 1.03, 1.41, 1.83 and 2.29. The values of f3 are slightly higher for the models with 
M*^*^ = 1.0 X 10^ and slightly lower for models with M<^^ = 2.0 x 10^ (Tabled]). The fraction of 
the enclosed mass of a merger object and the initial CC mass, ^ucf) is plotted as a function of /? 
in Figure [9l The merger objects of the models for the three initial CC masses show the same linear 
correlation with /3. 

The models with R^^ = 75 pc start already at the peri-galactic distance with a /3 of about one. 
Consequently, all 20 clusters merge. For the more extended models a number of star cluster are 
initially located outside the tidal radius. As we use a highly eccentric orbit, the CCs move rapidly 
towards larger Galactic distances, e.g. after 100 Myr they are at a distance of about 35 kpc. Due 
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to the lower gravitational field of the Milky Way at larger distances, the tidal radii increase leading 
to lower values of (3. Figure [10] shows the variation of fi with time for models with M^^ = 1.5 
X 10^. For models M_l_1.5_100, M_l_1.5_125, and M_l_1.5_150 the period with /3 > 1 is getting 
longer, resulting in lower numbers of merged clusters for larger CCs. The turnover in r^^, as shown 
in Figures Wp and El occurs at those R'^y^ i where /3 is sufficiently large to allow entire star clusters 
to escape the merging process. 

3.3. Variation of the initial distribution of star clusters 

We use the CC parameters of model M_l_1.5_100 as a basis to analyze the influence of the 
detailed distribution of star clusters in the CC. We calculate the evolution of five additional models, 
of which two (M_2_1.5_100 and M_3_1.5_100) have a similar concentration of clusters in their center 
as M_l_1.5_100, whereas the other three models (M_4_1.5_100, M_5_1.5_100 and M_6_1.5_100) show 
a less concentrated distribution of star clusters (see Figure Ht) . All six models result from exactly 
the same Plummer model but with different random number seeds. 

The average and the standard deviation of the effective radii and enclosed masses of the merger 
objects of all six models are r^^ = 23.8±3.2 pc and Menci = 0.95±0.17 x 10^ M©, respectively. 
The standard deviations, which correspond to relative deviations of 13% (refj) and 18% (Mend), 
provide an order of magnitude estimate of the influence of the distribution of star clusters on the 
structural parameters of the merger objects. 

The merger objects resulting from compact initial configurations 1 to 3 have on average a 
significantly higher mass (Mend = 1.10ib0.06 x 10^ Mq) than the less concentrated configurations 
4 to 6 (Mend = O.BlibO.OG X 10^ M0). In contrast, no clear difference in the effective radii of the 
merger objects can be seen between concentrated and extended initial distributions. 

For comparison, an additional model M_7_1.5_200 has been calculated. It has a CC Plummer 
radius of 200 pc and a relatively small cutoff radius of 400 pc. A scaled version of model M_l_1.5_100 
with R^y" — 200 pc and a CC cutoff radius of 800 pc, which is four times the CC Plummer radius, 
would have 6 star clusters beyond 400 pc. Model M_7_1.5_200 has the same cutoff radius and 
therefore the same values of (3 as model M_l_1.5_100. Due to the broader initial distribution of star 
cluster in M_7_1.5_200 the enclosed mass and the effective radius are slightly lower than in model 
M_l_1.5_100. 

The results demonstrate that next to the CC mass and the CC Plummer radius also the cutoff 
radius and the exact distribution of star clusters in a CC are key parameters for extended models, 
considerably influencing the structural parameters of the merger objects. 
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3.4. Comparison with Observations 



The enclosed masses were converted into absolute ^-magnitudes to allow for direct comparison 
with the observed data, using the formula 



My = Mv,solar " 2.5 X log(Mencl 



If' 



(7) 



where, Mv,soiar = 4.83 mag is the absolute solar V magnitude, Monci the enclosed mass of the merger 
obj ect and the mass-to -light ratio. We use a mass-to-light ratio of 2.05ib0.50, as determined 



by iBaumgardt et al.l ([20091) for NGC 2419. Apparen t V-magnitudes ca n by calculated by adding 
the distance modulus of NGC 2419 (m-M)o = 19.6 (IRipepi et al.l 120071 ^ and a V-band extinction 
of Ay = 0.25 using a reddening of E(B-V) = 0.08 ijRmem et ahlboOTl ). 



Figure [TT] shows an exemplary contour plot of the merger object of model M_l_1.5_100 projected 
onto the sky using Galactic coordinates. The lowest contour line corresponds to 32 mag arcsec"^. 
At this low surface brightness, NGC 2419 is detectable up to radii of about 17' corresponding to 
415 pc. Observed photometry, however, covers only a region within the inn er 60'' of NGC 2419 , 
while surface densities need to be derived from star counts in the outer parts. iRipepi et al.l (|2007l ) 
found evidence for stars associated with NGC 2419 up to radii of 15'. 

Figure [12] shows the effective radius, r^s, versus the absolute magnitude. My, of the observa- 
tions and the models as given in Table [1] and [3l respectively. Due to the scatter of the observed 
parameters of NGC 2419 and the turnover in the effective radii of the models as described in Section 
13. 2[ a large number of models are compatible with the observed parameters of NGC 2419. Only 
extended models with initial CC masses of M^*^ = 1.0 x 10® and compact models with M*^^ = 
2.0 X 10® are clearly incompatible with observations. 

A more detaile d comparison bet ween models and observations can be achieved using surface 
brightness profiles. [Bellazzinj ([2007[ ) compiled a surface brightness profile out to projected radii 
of 477'.'5. However, as he used radii between 720" and 900", which are well inside NGC 2419, 
as a reference field to estimate the background level, Bellazini's surface brightness estimates are 
expected to be slightly too faint in the outer parts NGC 2419. 

The radial surface brightness profiles of four exemplary models projected onto the sky and 
converted to units of mag arcsec"^ are shown in Figure [T3[ The merger objec ts show a King-like 
profile out to radii of about 1000". The observed profile from [Bellazzinj (|2007[ ) is added to Figure 
[T3] to allow for a direct comparison. The surface brightness profile of model M_l_1.0_50 agrees 
very well with the observed profile at all radii. The other three models in Figure [13] illustrate how 
the profiles change when one of the parameter mass, size and initial distribution of star clusters is 
modified. Model M_l_1.5_50 has a similar shape as M_l_1.0_50, but is too bright at all radii due to 
the larger mass. Model M_l_1.5_100, which is a more extended version of model M_l_1.5_50, agrees 
well with the observations between 10"and 100", but it is considerably brighter in the center and 
the outer parts. Model M_4_1.5_100, which has a broader initial distribution of star clusters than 
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model M_l_1.5_100, shows a smaller deviation from the observed surface brightness profile than 
model M_l_1.5_100. 



Baumgardt et al.l (|2009l ) observed the radial velocities of 40 stars within a projected radius of 
100 pc of NGC2419 and derived a velocity dispersion of o" = 4.14zt0.48 km s~^. The line-of-sight 
velocity dispersions within a projected radius of 100 pc of the models are listed in Table [3l A 
number of models with masses M*^^ = 1.0, 1.5, and 2.0 x 10^ Mq have velocity dispersions that 
are within the one sigma error of the observed velocity dispersion. Model M_l_1. 0_50, which has an 



effecti ve radius and enclosed mass very close to the newest observed values from iBaumgardt et al 



(|2009l ) and a surface brightness profile that is a good approximation of the observed profile, has a 



velocity dispersion of a = 4.13 km s ^, i.e. almost exactly the observed one. 

Considering masses, effective radii, surface brightness profiles, and velocity dispersions, model 
M_l_1.0_50 provides the best representation of NGC 2419. However, a number of models reproduce 
the observed structural parameters of NGC 2419 quite well within the observational uncertainties, 
demonstrating that an object like NGC 2419 can be formed from merged CCs without the need of 
fine-tuning of the input parameters. 



4. DISCUSSION 

4.1. The Merged Cluster Complex Scenario 

The proposed formation scenario for NGC 2419 starts with newly born complexes of star clus- 
ters in the Galactic halo with orbital parameters allowing for an highly eccentric orbit. We model 
the dynamical evolution of various CCs leading to merger objects. We do not, however, consider 
the process which formed the CCs in the first place as this would increase the complexity of the 
simulations and add more degrees of freedom making the interpretation of the results difficult. 

For all 27 models the majority of star clusters merge into a stable object. The turnover in the 
reff vs. Menci space (Figure [Tj M and [T2]) leads to degenerate states, because a relatively compact 
CC can produce a comparable merger object as a more massive CC having a significantly larger 
CC size. In consequence, a range of initial conditions can form a merger object comparable to 
NGC 2419 preventing us to pinpoint the parameters of the original CC, which formed NGC 2419. 
On the other hand, the larger the range of initial conditions that end up in a NGC 2419-like object, 
the larger is the probability of creating a massive EC like NGC 2419. 

As the individual star clusters of the CC formed at approximately th e same time from in olecular 



clouds of a galaxy, the observed absence of multiple stellar populations (iRipepi et al.ll2007l ) is fully 
consistent with our model. A small scatter in metallicity would also be explained, as the individual 
pre-cluster cloud cores of the complex could have had slightly different metallicities. In addition, 
some stars from the parent galaxy might have been c aptured by NGC 2419 during its formation as 



demonstrated for massive star clusters like u Cen bv lFellhauer et al.l (|2006l ) 
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Observations demonstrate that massive complexes of star clusters nowadays predominantly 
form during gravitational encounters between late-type galaxies. Young massive CCs with masses 
above 10^ Mq and sizes of a few hundred pc, comparable to those used as initial conditions for 
the numerical simu lations presented i n Sect. [3l have b een observed in the Antennae galaxies 
(|Bastian et ahlbood ) and in NGC922 dPellerin et al.lboiol '). 



Since galaxy-galaxy mergers are anticipated to have been much more common during early 
hierarchical structure formation it is expected that star formation in cluster complexes has been a 
significant star formation mode during early cosmological epochs. 



4.2. Potential Association with Stellar Streams in the Milky Way halo 



As stellar structures in the outer halo are long lived features, remnants of a parent galaxy or 
stellar tidal streams may still be observable in the Milky Way halo. 



Newberg et al.l (120031 ) found an over-density of A-type stars at a distance of 83 to 85 kpc, 
which has a width of at least 10° and which was traced for more th an 20° on the sky. N GC 2419 
is located within this debris (on the sky and at the same distance). I Newberg et al.l (j2003l ) argued 
that NGC 2419 and the stellar over-density might be associated with the Sagittarius dwarf galaxy, 
as NGC 2419 lies only 13 kpc from its orbital plane. However, the stellar streams fror n Sgr dwarf 
have mean metallicities between [Fe/H] = -0.4 near the core and -1.1 within the a rms (jChou et al 
200?! ). while NGC 2419 has a very low metallicity of [Fe/H] = -2.12 (|Harrislll99d ). 



The Vir go Stellar Stream (VSS) has a metallicity of [Fe/H] = -1.86 with a large scatter 
of 0.40 dex (IDuffau et al.l l2006l ). which makes it comparable to the metallicity of NGC 2419. 
Casetti-Dinescu et al.l (|2009l ) determined the proper motion of the VSS and calculated an orbit 
with a pericentric distance of 11 kpc and an apocentric distance of 89 kpc. They concluded that 
the current position of NGC 2419 is compatible with this orbit. 



Newberg et al.l (|2009l ) found a new stellar stream at a distance of 35 kpc in the constellation 



Cetus having a nearly polar orbit. The so-called Cetus Polar Stream (CPS) has a very low metal- 
licity of [Fe/H] = -2.1. While the metalhcity of CPS is in excellent agreement with NGC 2419, the 
polar orbit along Galactic longitude I = 143° appears to be in conflict with the current location of 
NGC 2419, which is at / = 180°. 



The Orphan Stream, inde pendently discovered b v lGrillmaid (120061 ) and lBelokurov et al.l ([200a), 
has been analyzed in detail bv lNewberg et al.l (j2010l ). They found a very low metallicity of [Fe/H] 
= -2.1 and a likely orbit between a pericentric distance of 16.4 kpc and an apocentric distance of 
about 90 kpc. These parameters are very similar to those of NGC 2419. 



Any potential association with the Sgr stream, VSS, CPS, or Orphan stream can only be 
verified after the proper motion of NGC 2419 has been measured. 
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4.3. NGC 2419-like Objects in other Galaxies 



NGC 2419 is located in the outer Galactic halo. Also the lower mass ECs of th e Milky Way are 
halo objects: 9 out of 13 ECs have galacto-centric distances greater than 20 kpc (jHarrislll996l ). A 
similar trend has been shown for the other two disk galaxies in the Local Group: 12 out of 13 ECs 
associated with M31 and both ECs found in M33 haye prgjected radii well o utside the optical disks 
of these galaxies (jHuxor et alJl2008l : IStonkute et alJl2008l : iHuxor et alJl2009l ). No EC as massive as 
NGC 2419 has been found in these galaxies. 

Due to the limited field of view of the Hubble Space Telescope, most extragalactic studies 
on GCs cover only (a part of) the optical disk of the respective galaxies. An outer halo object 
like NGC 2419 would be found only by chance if it is projected onto the main body of a galaxy. 
Another reason for incompleteness is the difficulty of distinguishing ECs from background galaxies. 
Therefore, a number of surveys applied a size limit to redu ce the contaminati on of background 
galaxies, e.g. the large GC surveys of the Virgo Cluster (jjordan et al.l l2005l ) and the Fornax 
Cluster (iMasters et al.l 120101 ) . covering dozens of galaxies. Their size limit of reff < 10 pc excludes 
all ECs from their GC catalogues. 

While NGC 2419-like objects with masses of Mec ^ 10^ Mq were not in the focus of extra- 
galactic surveys, much effort has been made to detect and to analyze stellar objects of M ^ lO'^ Mf?) 



since the di scovery of ultra-co i npact dwarf galaxies (UCDs) in the Fornax Cluster by lHilker et al 
(| 19991 ) and iDrinkwater et al.l (|2000l ) . A number of massive ste llar objects with e f fective radii of 
rpff K . 20 pc have been found i n the Fornax and the Virgo cluster (jMieske et al.ll2008l : iHasegan et al 



20051; lEvstigneeva et al.ll2007l V 



Richtler et al.l (|2005l ) present a NGC 2419-like object in the halo of the ehi ptical galaxy NGC 13 99 



which is the central galaxy of the Fornax Cluster. This object, labeled 90:12 in lRichtler et al.l (j2005l ). 
is located at a projected distance of about 40 kpc and has a very blue color indicating a very low 
metallicity. It has an absolute V-band magnitude of My = —10.04 mag, which corresponds to a 
mass of M90:i2 = 1.8 x 10^ M©, assuming the same mass-to-light ratio as NGC 2419. The effective 
radius of r^s = 27 pc is larger than that of NGC 2419. In Figure [121 this EC is located between 
the values of models M_l_2.0_100 and M_l_3.0_100. Therefore, the parameters of EC 90:12 are 
consistent with being a more massive version of NGC 2419. 

Considerably larger and more complete datasets of stellar halo objects like ECs and UCDs are 
necessary to draw statistically significant conclusions on the question whether NGC 2419 is a low- 
mass UCD, or a high-mass EC, and whether they are all remnants of merged CCs. A parametric 
study covering the entir e mass range from ECs to UCDs (M*^^ = 10^-^ to 10^ M©) wih be presented 



in a subsequent paper (IBriins et al.l 12011 
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5. SUMMARY 

The Galactic globular cluster NGC 2419 has unique parameters. It is one of the most luminous, 
one of the most distant, and as well one of the most extended GCs of the Milky Way. Apart from 
these unusual parameters, NGC 2419 appears to be a normal Galactic GC having a low metallicity 
and a single stellar population. 

We propose a new formation scenario for NGC 2419, being a remnant of a merged star cluster 
complex, which was formed during an interaction between a gas-rich galaxy and the Milky Way. 
To test this hypothesis, we performed particle-mesh-code computations of 27 different CC models. 
We vary the CC mass, the CC size, and the initial distribution of star clusters in the CC to analyze 
the influence of these parameters on the resulting obj ects. These CCs are comparable to those 



observed in the Antennae galaxies (iBastian et al.ll2006l ) and in NGC 922 (jPellerin et al.ll201Cl ). 



For all 27 models, the vast majority of star clusters merged into a stable object. The derived 
parameters mass, absolute V-band magnitude, effective radius, velocity dispersion and the surface 
brightness profile are, for a number of models, in good agreement with those observed for NGC 2419. 

The effective radii of the merger objects increase with increasing CC size up to = 75 pc 
for M^'~' = 1.0 X 10^ and up to 100 pc for the more massive models. For larger values of the 
effective radii decrease rapidly (Figure [7[). The turnover in the reg vs. Mgnci space (Figures [8] and 
[T2]l leads to degenerate states, as relatively compact CCs can produce a comparable merger object 
as a more massive CC having a significantly larger size. Despite the large range of CC sizes {Rpi" 
= 25 to 300 pc) of the models with initial CC masses of M'^'^ = 2.0 X 10^ the effective radii of 
the merger objects are constrained to values between 13.2 and 25.6 pc. 

In consequence, a range of initial conditions can form a merger object comparable to NGC 2419 
preventing us to pinpoint the parameters of the original CC, which formed NGC 2419. On the other 
hand, the larger the range of initial conditions that end up in a NGC 2419-like object, the larger is 
the probability of creating a massive EC like NGC 2419. 

Due to the limited field of view of the Hubble Space Telescope, most extragalactic studies 
on GCs cover only (a part of) the optical disk of the respective galaxies. An outer halo object 
like NGC 2419 would be found only by chance if it is projected onto the main body of a galaxy. 
Another reason for incompleteness is the difficulty of distinguishing ECs from background galaxies. 
Therefore, a number of surveys applied a size limit of about 10 pc to reduce the contamination 
of background galaxies. Thereby these surveys exclude all NGC 2419-like ECs from their GC 
catalogues. As massive galaxy-galaxy interactions are expected to have been more numerous in the 
past, a large number of NGC 2419-like objects probably await their detection in the outer halos of 
various galaxies. 

We conclude that NGC 2419 can be well explained by the merged cluster complex scenario. 
Measurements of the proper motion of NGC 2419 are indispensable to further study the proposed 
scenario and to potentially associate NGC 2419 with one of the stellar streams in the outer Galactic 
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Fig. 1.— Overview of GCs of the Milky Way (|Ha iris| [1999) . the LMC (jMackey fc Gilmorel 



2004 



van den Bergh k MackevI I2004I'). M51. and M81 (|Chandar et alJ 12004 ) . NGC 5128 (IHarris et al 



2OO2I). 68 dwarf g alaxies (IGeordev et all l2009l ) and t he ECs from M31 (IHuxor et al.1 12004 : 



Mackev et alJl2006l). M33 (IStonkute et al.ll2008l ). Scl-dEl (|Da Costa et alJl2009l l. and one EC from 



NGC 1399 ([Richtler et alj|2005l ) in the r^Q vs. My space. NGC 2419 is in a rather isolated position. 
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Fig. 2. — The orbit of NGC2419 traced back from its current position and projected to the xy-, 
the XZ-, and the yz-plane The stars indicate the observed current position of NGC 2419. Lower left: 
The distance of NGC 2419 to the Galactic center. The orbital period is about 1.3 Gyr. 
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Fig. 3. — The parameter range covered in the M^*-" vs. -Rpj^ space. The arrows indicate the increase 
of the parameters a and /? (see Sect. 13. 2p . The circle with the additional asterix marks the values 
of M^^ and of those models, where the initial distribution of star clusters within the CC was 
varied (see Sect. 13. 3p . 
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Fig. 4. — Three exemplary initial distributions of star clusters (small circles with radius i?cut) '^^ 
a CC (surrounding circle with radius Rc^) projected onto the x-y-plane. Model M_l_1.5_50 (b) 
is a scaled version of M_l_1.5_100 (a) only differing in Plummer radius R^i"- Some star clusters 
do already overlap in the center at the beginning. Model M_4_1.5_100 (c) has a less concentrated 
distribution of clusters than M_l_1.5_100. 



- 23 - 




Fig. 5. — Time evolution of the merger object in model M_l_1.5_100. Contour plots on the x-y- 
plane displayed at t' = 0, 50, 100, 300, 760 and 1500 Myr. The lowest contour level corresponds to 
5 particles per pixel. The pixel size is 5 pc. This yields 0.15 pc~^. The contour levels increase 
further by a factor of 3. 
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Fig. 6. — Left: Enclosed mass versus time after 1.5 Gyr of evolution of model M_l_1.5_100. The 
curve becomes fairly flat with increasing time. An almost stable merger object forms suffering only 
slightly from mass loss. Right: effective radius against time after 1.5 Gyr of evolution of model 
M_l_1.5_100. The effective radius decreases with time up to about 7 Gyr. Thereafter it is almost 
constant. 
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Fig. 7. — a: Enclosed mass, Mend, of the merger objects as a function of the initial CC Plummer 
radius, Rpi"- b: Effective radius, r^s, of the merger objects as a function of the initial CC Plummer 
radius, R'^p. 
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Fig. 8. — Effective radius r^s as a function of the enclosed mass of the merger objects for the CC 
models with R'^f^ = 25, 50, 75, 100, 125, and 150 pc (symbols from right to left) and M^^ = 1.0, 
1.5, and 2.0 x 10*^ M©. Models with the same initial CC mass are connected by solid lines. Models 
with the same R^^" are connected by dotted lines to illustrate the dependence on the initial CC 
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50, 75, 100, 125, and 150 pc for the first 200 Myr (the orbital period is about 1.3 Gyr). 
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Fig. 11. — Contour plot projected onto the sky in Galactic coordinates for model M_l_1.5_100 at 
the current observed position of NGC 2419. The pixel size is 0'.6. The contour levels go from 32 to 
20 mag arcsec"^ in steps of two mag arcsec"^. 
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Fig. 12. — Effective radius, r^s, as a function of absolute magnitude, My, for all computations (see 
Table El the dotted line connects models with R^y" = 100 pc) and observed values of NGC 2419 
(black stars). The error bars for the o bserved values are fr om the respective papers (see Table [1]). 
In addition, EC 90:12 from NGC 1399 (|Richtler et ahiboosl l is included (white star). 
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Fig. 13. — Surface brightness profiles of four exemplary models at the current position of NGC 2419 
using a mass-to-light ratio of 2.05. A V-band extinction of Ay = 0.25 was ap plied to th e mod els 
to allow for a direct comparison with the observed surface brightness profile of iBellazzinil (120071 ). 



- 32 - 



Table 1. Observational Parameters of NGC 2419 





BC^ 


Bb 


MvdM'^ 




Absolute V magnitude My [mag] 


-9.28 


-9.4 


-9.31 


-9.58 


Total mass'' [10^ M©] 


0.90 


1.01 


0.93 


1.19 


Central surface brightness /Uy [mag arcsec"^] 


19.61 


19.55 


19.44 


19.83 


Core radius rc [pc] 


7.8 


7.8 


7.9 


8.57 


Effective radius r^s [pc] 


19.2 


23.5 


19.9 


17.88 


Tidal radius rt [pc] 


190 


174 


204 


214 



Baumgardt et al.1 (|2009l ) 



Bellazzinil (I2OO7I ) 



McLaughlin van der Marell (|2005l ) 



fearrisi 



''using the mass-to-light ratio of 2.05 from iBaumgardt et al.l (|2009l ) 



Table 2. Initial Cluster Complex and Star Cluster Parameters 



Parameter range 

Star Cluster (SC) 

Number of star cluster particles, A'g 
Initial SC mass, M^c [iqS M©] 
Plummer radius of the SC, i?pp [pc] 
Cutoff radius of the SC, [pc] 
Cluster Complex (CC) 
Number of star clusters, Nq'^ 
Initial CC mass, M*^^ [10^ M©] 
Plummer radius of the CC, R^^^ [pc] 
Cutoff radius of the CC, R^-^ [pc] 
a-parameter 



100 000 
0.5 - 1.5 
4 
20 



20 

1.0 - 3.0 
25 - 300 
100 - 1200 
0.160 - 0.013 
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Table 3. Parameters of the Merger Objects 



ModeP 


' iVi 


-^^^^ciicl 




P 

' 11 


f 

' CII 


Tn ^ 
' C 




a' 






(10*^ Mn) 


(maEl 


(pc) 


(pc) 


(pc) 


(mag arcsec"^) 


(km s~^) 


M_l_1.0_25 


20 


0.93 


-9.31 


18.8 


14.2 


6.6 


18.99 


4.88 


M_l_1.0_50 


20 


0.85 


-9.21 


25.8 


19.2 


5.7 


19.36 


4.13 


M_l_1.0_75 


20 


0.77 


-9.11 


29.8 


21.8 


5.1 


19.41 


3.80 


M_l_1.0_100 


19 


0.64 


-8.91 


24.9 


19.1 


5.5 


19.51 


3.59 


M_l_1.0_125 


16 


0.54 


-8.72 


22.1 


16.6 


5.0 


19.47 


3.51 


M_l_1.0_150 


14 


0.46 


-8.55 


19.6 


14.4 


4.5 


19.39 


3.47 


M_l_1.5_25 


20 


1.38 


-9.74 


19.9 


14.9 


6.7 


18.65 


5.85 


M_l_1.5_50 


20 


1.29 


-9.67 


26.6 


20.0 


6.2 


19.04 


4.99 


M_l_1.5_75 


20 


1.17 


-9.56 


30.1 


22.5 


5.6 


19.12 


4.57 


M_l_1.5_100 


19 


1.05 


-9.44 


31.8 


23.5 


4.9 


19.08 


4.34 


M_l_1.5_125 


17 


0.92 


-9.30 


28.8 


21.3 


5.3 


19.15 


4.21 


M_l_1.5_150 


15 


0.80 


-9.15 


25.1 


18.8 


4.8 


19.09 


4.10 


M_l_2.0_25 


20 


1.84 


-10.05 


21.1 


15.9 


7.2 


18.47 


6.55 


M_l_2.0_50 


20 


1.75 


-10.00 


28.4 


21.2 


6.4 


18.83 


5.65 


M_l_2.0_75 


20 


1.61 


-9.91 


32.0 


23.7 


6.0 


18.88 


5.26 


M_l_2.0_100 


19 


1.48 


-9.82 


35.2 


25.6 


5.3 


18.87 


4.99 


M_l_2.0_125 


16 


1.26 


-9.64 


31.5 


23.2 


5.4 


19.01 


4.76 


M_l_2.0_150 


18 


1.17 


-9.56 


28.9 


21.2 


5.1 


18.81 


4.81 


M_l_2.0_200 


14 


0.99 


-9.38 


26.4 


18.6 


4.7 


18.85 


4.64 


M_l_2.0_300 


10 


0.62 


-8.87 


17.6 


13.2 


4.8 


19.03 


4.18 


M_l_3.0_100 


19 


2.40 


-10.34 


39.6 


28.8 


5.3 


18.45 


6.08 


M_2_1.5_100 


19 


1.17 


-9.56 


37.2 


27.5 


5.8 


19.34 


4.24 


M_3_1.5_100 


19 


1.07 


-9.46 


31.9 


23.7 


5.3 


19.21 


4.33 


M_4_1.5_100 


18 


0.80 


-9.15 


22.7 


17.9 


5.8 


19.23 


4.10 


M_5_1.5_100 


19 


0.87 


-9.24 


31.3 


24.5 


5.4 


19.35 


3.84 


M_6_1.5_100 


15 


0.76 


-9.09 


33.4 


25.5 


5.2 


19.58 


3.53 


M_7_1.5_200 


19 


0.93 


-9.31 


30.1 


21.9 


5.1 


19.12 


4.22 



''Model__Configuration__M*"'"__i?^ 
^Number of merged star clusters. 
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'^Enclosed mass of merger object within 800 pc. 

Absolute V-magnitude of merger object. 
^Half-mass radius of merger object. 

^Effective radius, i.e. tlic projected half-mass radius of the merger object. 
sCore radius of merger object obtained from a King fit. 

^Central surface brightness fiy of merger object obtained from a King fit. An extinction 
of Ay = 0.25 has been applied. 

'Velocity dispersion within a projected radius of 100 pc. 
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Tablc 4. /3-vahios of oiir initial CC models 





25 


50 pc 


75 pc 


100 pc 


125 pc 


150 pc 


200 pc 


300 pc 


1.0 X 10^ Mq^ 


0.38 


0.77 


1.19 


1.63 


2.14 


2.73 






1.5 X 10^ M0 


0.33 


0.68 


1.03 


1.41 


1.83 


2.29 






2.0 X 10^ Mq 


0.31 


0.62 


0.94 


1.28 


1.64 


2.04 


3.02 


7.27 


3.0 X 10^ Mq 








1.11 











^Plummer radius R^^ of the CC. 
^Initial CC mass 



